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1.  Introduction 

While  liquid  solution  reactions  are  much  more  important  in  chemistry,  gas  phase  reactions  are 
much  better  understood.  Given  the  central  importance  of  solution  reactions  to  inorganic,  organic, 
industrial  and  biochemistry,  it  is  rather  surprising  that,  as  yet,  there  is  not  a  single  such  reaction 
whose  molecular  dynamics  are  understood  in  detail.  Theoretical  and  experimental  evidence  already 
makes  clear  that  much  of  the  important  molecular  dynamic  action  in  solution  reactions  occurs  on 
the  picosecond  and  subpicosecond  time  scales.  The  dihalogen  photodissociation  and  recombination 
reactions,  Xi  +  hv  —  X  4-  X— X2 ,  involving  the  simplest  possible  molecular  reactants  and  pro¬ 
ducts,  diatomics,  and  in  rare  gas  solution  involving  only  two  elements,  seem  excellent  candidates 
for  study. 

2.  Theory 

The  first  deterministic  theoretical  study  of  the  molecular  dynamics  of  reactions  was  by  BUNKER 
and  JACOBSON [1],  who  computed  the  classical  trajectories  for  I2  in  CC14  solvent  represented  by  26 
spherical,  structureless  particles  in  a  specular  cube.  MURREL,  ST  ACE  and  DAMMEL[2] 
modelled  the  photodissociation  of  I2  in  dense  inert  gases,  I2  plus  22  gas  atoms  in  a  spherical,  soft- 
walled  container.  We  have  similarly  modelled  I2  plus  50  Xe  atoms  at  liquid  density  in  truncated 
octahedral  periodic  boundary  conditions [3, 4),  computing  the  photodissociation,  solvent  caging, 
atomic  recombination  and  vibrational  energy  decay  to  the  solvent  from  the  new  I2  molecule,  as 
shown  in  Fig.  1.  The  conclusion  of  all  three  molecular  dynamic  studies  is  that  geminate  recombina¬ 
tion  is  usually  a  very  fast  process,  over  within  a  few  picoseconds.  An  important  caveat,  and  a 
weakness  in  these  theoretical  studies,  is  that  the  process  whereby  the  I  atoms  dissociating  on  an 
excited  state  potential  surface  refind  the  ground  state  surface  on  which  they  recombine  is  not  well 
understood,  and  is  therefore  handled  in  these  calculations  by  arbitrary  assumptions  which  may  be 
incorrect.  If  so,  the  real  time  for  geminate  recombination  may  be  longer  than  the  few  picoseconds 
calculated. 

Fig.  1.  Time  evolution  in  liquid  Xe  solution  of 
I2  vibrational  energy  during  reaction  sequence 
of  photodissotiation  (at  time  zero),  solvent  cag¬ 
ing  of  some  of  the  recoiling  I  atoms,  radical 
recombination,  and  vibrational  loss  to  solvent. 
Also  shown  is  the  equilibrium  Boltzmann  vibra¬ 
tional  energy  distribution  before  photodissocia¬ 
tion.  The  vibrational  enetgy  distribution  al¬ 
ready  bifurcates  into  two  branches  within  the 
0-5  ps  period,  the  higher  one  corresponding  to 
those  I  atom  pairs  which  have  escaped  the  cage 
and  whose  minimum  energy  is  the  I2  dissocia¬ 
tion  energy,  and  the  lower  one  corresponding 
to  recombined  I2  progressively  losing  vibration¬ 
al  energy  to  the  solvent. 
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In  addition,  theoretical  calculations  by  NESBITT  and  HYNES(S,6]  for  I2  in  rare  gases  and  in 
CQ4  and  by  our  group (3, 4]  in  liquid  Xe  indicate  that  the  decay  to  the  solvent  of  the  vibrational 
energy  in  the  newly  reformed  I2  molecule  will  require  the  order  of  a  hundred  to  hundreds  of 
picoseconds,  as  shown  in  Fig.  1,  considerably  slower  than  the  time  required  for  geminate  recombi¬ 
nation.  Note  that  these  calculations  are  for  solvent  atoms  or  molecules  which  are  very  weakly 
bound  to  one  another,  and  that  the  vibrational  relaxation  might  be  quite  different,  for  example,  in  a 
strongly  hydrogen-bonded  liquid. 

From  our  calculated  molecular  dynamics,  plus  the  potential  curves  and  transition  dipole 
moments,  we  can  compute  transient  electronic  absorption  spectra[7],  as  shown  in  Fig.  2,  which 
includes  the  A'—X  ,  B—X,  and  the  B'  luf'n)  —  X  transitions.  A  related  nonmolecuiar  dynamics 
spectral  calculation  has  been  carried  out  by  NESBITT  and  HYNES [6].  A  small  quantum  correction 
by  temperature  scaling,  which  would  be  exact  for  the  coordinate  distribution  in  the  harmonic  limit, 
is  made  to  the  equilibrium  spectra  which  agree  well  with  the  known  I2  experimental  spectral 
points(8].  Note  that  spectra  measured  at  different  wavelengths  follow  different  time  histories, 
which  can  in  principle  be  used  to  follow  the  time  evolution  of  the  vibrational  energy  distribution  of 
the  relaxing  I2  n  olecules. 


WAVELENGTH  (nm) 


ENERGY  (loW1) 


Fig.  2.  Transient  electronic  absorption  spectra 
computed  from  molecular  dynamics  for  I2  (in 
Xe  solution)  reaction  sequence  of  photodissoci¬ 
ation,  solvent  caging,  radical  recombination, 
and  vibrational  decay.  Time  zero  is  the  photo¬ 
dissociation.  The  top  curve  is  the  computed 
equilibrium  spectrum  before  photodissociation 
on  which  are  superimposed  the  experimental 
points  from  Tel!inghuisen[8]. 


3.  Experiment 

The  earliest  picosecond  experimental  results  were  by  the  EISENTHAL  group(9, 10],  who  measured 
the  transient  electronic  absorption  spectra  after  excitation  at  530  nm.  Decay  times  of  ~70  ps  for  I2 
in  hexadecane  and  —140  ps  in  CCI4  were  observed.  Subsequently  these  studies  were  extended  by 
the  EISENTHAL  group  to  I2  in  aromatic  solvents]!!]  which  are  believed  to  form  complexes  with  1 
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atoms,  and  by  LANGHOFFI12]  who  observed  1}  photodissociation  in  several  weakly  associated 
liquids,  finding  decay  times  in  the  -100- 150  ps  range.  More  recently,  KELLEY  and 
RENTZEP1SU3)  have  observed  l2  photodissociation  in  fluid  and  liquid  Xe  with  a  decay  time  or 
—40  ps,  as  well  as  in  CCU,  and  similar  experiments  have  been  carried  out  by  the  PETERS 
group|14].  All  of  the  above  studies  used  second  harmonic  Nd  pump  light  at  —  530  nm.  which,  at 
least  in  the  gas  phase,  results  in  excitation  largely  to  the  bound  B  0u4(3fl)  state  which  presumably 
predissociates,  but  may  absorb  another  photon  in  the  meantimeUS,  16).  A  delay  of  —20  ps 
between  excitation  and  the  maximum  in  the  absorption  curve  has  been  attributed  to  absorption 
from  the  B  state  [9, 13, 17].  Thus,  the  linking  of  these  experimental  transient  spectra  to  the  molecu¬ 
lar  dynamics  of  the  I2  photodissociation  and  recombination  reaction  is  made  difficult  by  the  prob¬ 
able  presence  of  the  additional  processes  of  B  state  predissociation  and  absorption.  We  have 
recently  reported  [3]  transient  absorption  for  I2  excited  using  an  Ar+  -  dye  (DCM)  synchronously 
pumped  source  at  710  nm  into  the  dissociative  A  state,  thus  avoiding  the  problem  of  predissocia¬ 
tion.  In  order  to  achieve  the  sensitivity  needed  to  detect  the  very  weak(8]  I2  A  state  absorption,  we 
use  a  multiple  modulation  system[18]  based  on  the  discovery  by  HERITAGE[19]  and  LEVINE  and 
BETHEA120, 21]  that  the  noise  in  Ar*  synchronously  pumped  dye  lasers  falls  off  by  several  orders 
of  magnitude  in  going  from  the  audio  to  the  radio  frequency  region.  We  modulate  the  pump  and 
probe  beams  at  two  different  radio  frequencies  and  detect  at  the  difference  frequency,  using  inex¬ 
pensive  and  readily  available  radio  amateur  equipment(18].  In  addition,  we  audio  modulate  the 
pump  beam  and  synchronously  detect  at  that  frequency.  Decay  times  are  shown  in  Fig.  3.  We  sug¬ 
gest  the  hypothesis  that  vibrational  decay  might  be  expected  to  be  faster  for  the  more  strongly 
hydrogen-bonded  solvents  which  are  expected  to  have  a  greater  spectral  mode  density  in  the  range 
of  I2  vibrational  frequencies.  All  the  above  experimental  measurements  suffer  in  interpretation 
from  the  rather  weak  connection  between  molecular  dynamics  and  electronic  absorption,  which  is 
further  complicated  by  the  possible  presence  of  I2- solvent ,  1- solvent,  and  12  — 12  complexes. 


Fig.  3.  Experimental  transient  electronic 
absorption  spectra  for  I2  in  ethylene  glycol 
(0.2  molar,  1.3xlO-2  mole  fraction),  in 
ethyl  alcohol  (0.33  molar,  1.9x  10“ 2  mole 
fraction),  and  in  CC14  (0.0%  molar. 
9.3xlO~3  mole  fraction).  All  the  spectra 
are  for  perpendicular  orientation  of  linearly 
polarized  pump  and  probe  beams. 
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4.  Discussion  and  Conclusion 

In  summary,  two  hypotheses  have  been  advanced:  i)  that  geminate  recombination  for  I2  is  relatively 
slow  and  thus  accounts  for  the  observed  range  of  transient  absorption  decay  times  and  ii)  that  gem¬ 
inate  recombination  is  relatively  fast  and  that  vibrational  decay  times  to  reach  vibrational  levels  with 
higher  absorption  instead  account  for  the  observed  transient  absorption  decay  times.  At  the  present 
time,  the  available  theoretical  and  experimental  tools  have  not  been  sufficiently  powerful  to  cleanly 
disprove  either  or  both  of  these  hypotheses. 

Improvements  in  both  theoretical  and  experimental  tools  should  lead  to  more  stringent  tests. 
On  the  theoretical  side,  molecular  dynamics  and  spectral  calculations  can  certainly  be  extended  to  a 
variety  of  different  solvents,  and  to  different  pump  photon  energies  and  thus  different  I  atom  recoil 
energies.  In  addition,  different  assumptions  as  to  the  mechanism  for  relaxation  to  the  ground  state 
potential  curve  can  be  tried  out.  Transient  electronic  absorption  measurements  suffer  from  the 
intrinsic  limitation  of  a  weak  connection  to  molecular  dynamics.  In  contrast,  transient  Raman  spec¬ 
tra  can  directly  reveal,  for  example,  rotational  and  vibrational  periods  in  the  evolving  reactants  and 
products.  We  have  computed  such  transient  spectra  (3,4]  for  this  reaction  sequence.  The 


equivalent  resonance  Raman  spectra  can  also  be  calculated  and  they  also  will  reveal  much  of  the 
underlying  molecular  dynamics,  possibly  including  upper  state  recurrence  times[22].  Thus  Raman 
spectra  could  provide  a  definitive  test  of  the  above  hypotheses. 

It  is  surprising,  but  true,  that  for  no  chemical  reaction  in  solution,  not  even  for  one  as  sim¬ 
ple  and  as  well  studied  as  I2  photodissociation  and  recombination,  are  the  detailed  atomic  motions 
by  which  it  occurs  yet  known.  Even  such  basic  aspects  are  not  certain  as  the  order  of  magnitude  of 
the  time  required  for  caging  and  geminate  recombination,  and  whether  a  direct  deterministic  or  a 
stochastic  diffusions!  approach  to  geminate  recombination  is  most  appropriate. 

This  situation  may  soon  change,  as  there  is  now  a  four  order  of  magnitude  time  range,  —100 
fs  to  1  ns,  over  which  molecular  dynamic  calculations  and  short  light  pulse  experiments  can  over¬ 
lap.  Transient  infrared,  Raman  and  electronic  absorption  spectra  all  reflect  the  underlying  molecu¬ 
lar  dynamics  of  chemical  reactions  and  can  provide  an  interface  at  which  theory  and  experiment 
may  meet.  By  comparing  transient  spectra  computed  from  molecular  dynamics  with  the  equivalent 
measured  spectra,  one  can  hope  to  discover  the  microscopic  dynamics  by  which  many  chemical 
processes  occur. 

Thanks  for  the  support  which  has  made  this  work  possible  to  NSF  Chemistry,  ONR  Chemis¬ 
try,  NASA- Ames,  NIH  Division  of  Research  Resources,  and  Fonds  National  Suisse  for  fellowship 
support  to  P.  Bado. 
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